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Abstract. The distribution engendered by successive splitting of one point vortex are considered. The 
process of splitting a vortex in three using a reverse three-point vortex collapse course is analysed in great 
details and shown to be dissipative. A simple process of successive splitting is then defined and the resulting 
vorticity distribution and vortex populations are analysed. 
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1 Introduction 

Turbulent flows are often characterized by cascade dynam- 
ics that explain the transfer of quantities (energy, enstro- 
phy) between scales [Tj. In particular, in two-dimensional 
turbulent flows, a striking features is the presence of an 
inverse energy cascade, which leads to the emergence of co- 
herent vortices, dominating the flow dynamics [2I3I4I5I6] . 
together with the direct cascade of enstrophy. In order to 
tackle these problems, point vortices have been commonly 
used with some success: indeed they can approximate the 
inviscid dynamics of finite-sized vortices [71819] . as for in- 
stance in punctuated Hamiltonian models |6I10| . In these 
models, the advection of well-separated vortices is approx- 
imated via Hamiltonian point-vortex dynamics (thus in- 
viscid i.e. dissipationless); to account for the change in 
the vortex population toward smaller number of bigger 
vortices, dissipative merging processes are then included 
for vortices which have approached each other closer than 
a certain critical distance. In decaying two-dimensional 
turbulence, the merging process results in fact from the 
interaction of a few number of close vortices (T3] so that 
the understanding of low dimensional vortex dynamics is 
an essential ingredient of the whole picture [T3] . More gen- 
erally, it has been shown that point vortices can exhibit 
both the features of extremely high-dimensional as well as 
low-dimensional systems |llj . 

Since the pioneering work of Onsager on two dimen- 
sional turbulence [TJ] , a statistical approach to turbulence 
using point vortices has also been developed. These vor- 
tex systems display negative temperature, corresponding 
to states where same-sign vortices bound to form larger 
vortices [15114] . although special care in the definition of 
the thermodynamic limit has to be done[16j. In the same 



spirit, stationary flows resulting from point vortices can be 
obtained [17118119] and different kinetic theories can also 
be derived (see for instance [20] and references therein). 

Regarding the dynamical aspect, the relevance of point 
vortices, seen as exact solutions of the Euler equation, is 
still debated: indeed finite time singularities are present 
in the coupled dynamics of many point vortices. These 
singularities correspond to the collapse of three vortices 
[21122123] , and can been seen as the consequences of an ill 
posed problem. Such singularity arises in fact naturally in 
the Hamiltonian point vortices dynamics. Regardless the 
influence of the viscosity that would become dominant at 
short time before the singularity, it is important to notice 
that the dynamics toward the three vortices collapses is 
in fact reversible. Thus, it is tempting to consider the re- 
verse dynamics which would consist of vortices separation 
that has to be present in this Hamiltonian system! Such 
dynamics has been omitted until now since the number of 
point vortices was taken constant or only decreasing. 

Therefore, in this paper we take a different perspective 
on the existence of this finite time singularity, using it as 
a potential source of vorticity and vortices. One of the 
perspective is to offer the possibility of a statistical me- 
chanics approach with a varying number of vortices. At 
this stage, we remain within the Hamiltonian framework 
of point vortices and in particular we do not consider the 
regularization due to the viscosity for real fluid. In addi- 
tion, we focus here on the statistical distribution of vortic- 
ity rather than on the coupled dynamics of the system of 
vortices issued of this splitting. To be more precise, in this 
paper our main goal is to describe the statistical properties 
emerging from this genuine process of vortex splitting: we 
consider the distributions of vortices and vortex strength 
resulting from this simple mechanism of successive split- 
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ting according to the reverse collapse of one point vortex. 
We shall refer to this result as the offsprings of a point 
vortex. 

The paper is organized as follows, first we recall briefly 
the notion of point vortices, and how they naturally ap- 
pear as a solution of Euler equation. We then consider the 
process of splitting one vortex in three. We introduce the 
relevant parameters and analyse briefly the preliminary 
consequences. In particular, we observe that this splitting 
is inherently a dissipative process. Finally we define simple 
rules for successive splitting. Computing the offsprings of a 
point vortex with these rules we deduce general properties 
for the vortices distribution in the limit of high splitting 
processes. 



2 Basic equations 

Point vortices are singular solutions of some bidimensional 
physical systems described by a conservation equation of 
what we shall call a generalised vorticity Q given by 



on 

— + {Q,i>} = Q 



where {•,•} denotes the usual Poisson bracket, and ip is 
a stream function. The actual relation J? = F(ip) may 
depend on the considered physical system. For instance 
for the Euler equation it is simply given by Q — — V 2 ^. 
When Q = —V 2 ip + tp/ p 2 where t/j is, in this context, 
related to the electric potential (in suitable units) in a 
plasma, p s is the hybrid Larmor radius. Point vortices are 
defined by a vorticity distribution given by a superposition 
of Dirac functions, 



1 N 



*(*)) > 



where r = (x, y) is a vector in the plane of the flow, ki is 
the strength of vortex i (circulation), N is the total num- 
ber of vortices, and Ti(t) is the vortex position at time t. 
Using this expression of the vorticity and solving the Pois- 
son equation, in the Euler case, or the Helmholtz equation 
in the more general case, one obtains the current function 
associated to the point vortices. Thanks to the Helmholtz 
theorem, the motion of the vortices is determined by the 
value of the velocity field created by the other vortices 
at the position of the vortex. The point vortex motion is 
Hamiltonian and given by 



dH 

dxi 



dH 



,(i = l,'--,N) (3) 



where the Hamiltonian H is given by 

H = — ^kikjUQri - r 3 



(4) 



with for an unbounded plane U(x) — — log(x) in the Euler 
case, and U(x) = Kq(x) in the more general case of the 



plasma model (when p a — > oo the modified Besscl function 
Ko tends to the logarithm). The Hamiltonian ((4]) exhibits 
clearly that a system of point vortices is invariant by trans- 
lation and rotation, which implies both the conservation 
of the centre of vorticity and the total angular momentum, 
the motion of three vortices is integrable, while for four 
or more vortices Hamiltonian chaos come into play. When 
the distance between the vortices is smaller than the typ- 
ical interaction length (p s ) the behaviour of the two sys- 
tems is similar, while in the large distance limit, the Ko 
interaction decreases exponentially and the vortices are 
almost free. In the following, analytical computations will 
be made using the logarithmic interaction, which corre- 
sponds to the Euler flow, and we can expect the results to 
be qualitatively valid for the Bessel interaction as long as 
the vortices are not "too far" from each other. 

Regarding the Hamiltonian (TJJ , we remind that it does 
not represent the energy of the fluid (which is infinite al- 
ready with one point vortex), but rather corresponds to 
an energy of interaction between the vortices and is the 
one that is traditionnally used when making a statistical 
physics approach of point vortex systems. It is however 
important to recall there is actually an infinite "reserve" 
of energy in these systems. In what follows, we abusively 
refer to this interaction-energy, as energy. 

Let us now focus on a situation with only three vor- 
tices. In this restricted situation, it is easier to tackle the 
motion of the vortices by studying in fact their relative 
motion. Namely the three vortices form a triangle, and the 
relative motion describes the deformation of this triangle 
121124125126125] . The invariance by translation of ©, al- 
lows us a free choice of the origin of the plane, which we 
choose to be the centre of vorticity (when it exists). The 
other constants of motion written in a frame independent 
form become, 



H = [kik 2 lni? 3 + A;ifc 3 lnR 2 + k 3 k 2 lnR 1 ] 



K = k x k 2 R% 



kihR 2 2 



k 3 k 2 Ri , 



(5) 



where Ri = r/b|, with i ^ j ^ k. In fact it has been 
known for a long time that the motion of vortices can 
lead to singular solutions and finite time singularities, the 
most striking one occurring for a system of three vortices 
resulting in the collapse of the vortices in a finite time 
[21 22 23J or by time reversal, to an infinite expansion of 
the triangle formed by the vortices. The collapse or infinite 
expansion of the three point vortices are obtained when 
the following conditions are satisfied 



K = 



(6) 
(7) 



the harmonic mean of the vortex strengths ([7]) and the 
total angular momentum in its frame free form are 
both equal to zero [25126127] . 

In this paper we shall use this specific singularity and 
consider vorticity distribution arising from successive split- 
ting of a point vortex according to the collapse conditions. 
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Note that other type of singularities involving more vor- 
tices are effectively possible, however for the sake of sim- 
plicity, we restrict ourselves to reverse three-vortex col- 
lapse rules. In what follows we shall refer to the successive 
splitting process as the point vortex offsprings. 



3 Splitting of a Point Vortex 
3.1 Splitting rules 

In order to be consistent with the physical properties of 
vortex collapse, we successively divide vortices according 
to the collapses rules and keep the total vorticity constant. 
The splitting rules from one generation n of a vortex i of 
strength ki jTl to the next generation n + 1 read 



i=l 



3 

^ k 

i=l 



1 



z,n+l 



(8) 



(9) 



These equations are equivalent to 



i=l 
3 

V k 2 - k 



1 

i.n ' 



(10) 

(11) 



i=l 



the ki yTl+ i lie on the circle at the intersection of the sphere 
of radius \ki^ n \ defined by Eq. (fTTj) and the plane defined 
by Eq. (| 10[) . We therefore discuss the splitting in terms 
of the vector k„ + i = (&i >n+ i, fc 2 ,n+i, fe.n+i)- Tne tip of 
the vector lies on a circle, hence we parametrise it with an 
angle 6 n as: 



k n+ i = k n a + 




0„u + sin0 n v) ) , (12) 



where a = (l,l,l)/3, and for instance u = (1,— 1,0)/a/2 
and v = \/3aAu. In fact the choice of 9 n can be restricted 
to be picked within the segment [0, 2ir/3], since Eq. (|12l) 
exhibits the symmetry k 3}Tl (9 + 2/3tt) = ki iTl , k 3 ^ n (9 — 
2/3n) = k 2>n and thus all configurations up to a rela- 
belling can be obtained within this interval. 



3.2 Is the splitting always possible? 

We imagine that we are dealing with a system with many 
vortices (2n + 1 after n splittings) and consider the split- 
ting of one point vortex according to the rules ^ and ([7]) 
in this system. The total number of vortices changes, see 
Fig [1] for an illustration of the process. 
Regarding the energy we have 



SH n = H n+ i — H n — — — \n.A n , 

4-7T 



(13) 



= -— [kik 2 \nR 3 + k 1 k 3 lnR 2 + k 3 k 2 lnR 1 ,(14) 

where for instance the vortex k^ — ki + k 2 + k 3 was split 
in three. In order to be dynamically compatible, we first 
neglect the influence of the other vortices (which are con- 
sidered far enough to not interfere locally), but still we 
need to make sure that there is at least one triangle satis- 
fying the conditions and define the value of 5H n associated 
with the splitting. For that purpose after the splitting we 
name 2 and 3 the vortices whose strengths have the same 
sign with |jfea| < j^l, exponentiating Eq. (flU)) gives 



A r 



(15) 



then we divide by one noticing that R± 3 = 1, and 
obtain 



(16) 



Since the collapse is scale free (self-similar), we choose R\ 
as our length units, note r 2 = R 2 /R\, r 3 = R 3 /R\ and 
arrive at 

-k 2 /k 3 (17) 



The condition 



r 2 - ^ n ' 3 

becomes once rescaled 



hk 2 rl + hhrj + k 3 k 2 =0, (18) 
combining this last expression with ([7]) we have 

fc 2 (r 2 -l)+fc 3 (rf-l)=0, (19) 



and noting a = k 2 /k 3 (< 1), we finally obtain 



arf 1+a) + A 2 J k ^ =(l+a)r : 



2a 
3 



(20) 



A n being positive, aX 1+a — (1 + a)X a shows that r 3 £ 

]0 ; 1 + l/a[ and that < A 2 J klk2 < 1. The equality A n = 
1 being reached for r 3 — 1 which implies an equilateral 
triangle. This last configuration has to be excluded as it is 
dynamically a fix point, i.e the triangle does not expand 
or shrink hence can not be a starting point for a vortex 
splitting. So , since k\k 3 < 0, this means that A n > 1 and 
that 5H n — — A^hiA n < 0. The splitting of a vortex in 
three is thus a dissipative process! 

We now enforce that the solution is a triangle, using 
the rescaled variables this means that 



l = ri 



2r 3 r 2 cos (p . 



which implies (using Eq. (fT9 



COS If 



1(1 



a 



(21) 



(22) 
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Fig. 1. Simple 5th order vortex lineage. In this picture 
we see that after 8 splittings, 4 vortices are remaining on 
level 3. 
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Fig. 2. Probability density function (PDF) of the abso- 
lutes values of vortex strengths after respectively (from 
right to left) 10 3 , 5 I0 3 , 10 4 , 2 40 4 , 5 10 4 , 10 5 vortex splits. 
The PDF have been averaged over 64 trees, for the large 
values and up to 512 trees for the smallest ones. The shape 
appears as self-similar. 



The right hand side of Eq. (|2"2"|) is always positive, so 
— 7r/2 < cp < 7r/2. The minimum is obtained for r| = 
a(l + a)/(l + a 2 ) and equal to a 1 / 2 / (1 + a), which is al- 
ways smaller than one. Therefore, there is always a range 
of possibilities available for r% and the splitting is always 
possible. Note that for a given value of r% the two mirror- 
ing shapes of the triangle are possible, one giving rise to 
expansion (splitting) the other one to collapse. 



4 Vorticity distributions 

We are now interested in the vorticity distribution we 
obtain from such process. For this purpose we compute 
different trees originating from one vortex of strength k — 
1 as depicted in Fig. [T] The distributions are computed 
by successive vortex splitting. After each division, a vor- 
tex is chosen randomly among the global population and 
is split according to the reverse collapse rules, the division 
is the result of the uniformly random choice of an angle 9 n 
(see Eq. (fT2"j) '). Other possible rules made by assigning dif- 
ferent probabilities on vortices will be explored in future 
work. Results of the obtained distributions are depicted in 
Fig. [5] We shall notice that the chosen rules gives rise to 
a large spectrum of vorticities (see the logarithmic scale 
in FigJ5]). We notice as well that as the number of divi- 
sion increases the distribution spreads and the location of 
its maximum is slowly moving towards smaller values of 
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Fig. 3. Variation of energy versus number of created vor- 
tices. The curve has been averaged over 100 trees. Split- 
tings have been performed up to N — 65537 vortices start- 
ing from one vortex of strength ko = 1. One can notice 
a small logarithmic decay of energy versus the number of 
vortices, of the type SE s» — Alog 10 (iV) with A ps 5.9 10 -3 . 



the vorticity, moreover, the evolution of the distribution 
appears to have some kind of self-similar behaviour. In 
the same spirit the variation of energy as a function of 
generated vortices can be monitored. Results are shown 
in Fig. [3] and show a slow logarithmic decay of the total 
energy. 

In order to analyse this in more details we need to 
characterise the lineage (see Fig. [TJ after a given number 
of splittings. We will note the total number of splittings 
that occurred n. These successive divisions generate a tree 
(the phase space of the process) with 3 ra leaves at the 
extremities. Each division results in the choice of an angle 
9i. Now let us consider a particular "lineage" of order n, 
it gives rise to a family of N = 2n + 1 vortices. At each 
step of the division process (from n to n + 1) we choose 
any already existing vortex and split it with the rules ([121 . 
The trajectory in phase space corresponds to a connected 
graph with 2n + 1 leaves starting from the top of the tree 
and of total length 3n + 1. The number of possible graphs 
on the tree after n splittings is: (2n + l)!/(2 n n!). To move 
further on, and due to the large amount of possible graphs, 
we consider the global occupation Mi of the level i (see 
Fig. [T|) . We note Mi (n) the average number of leaves at 
level i at time n. Then we have 



, , Mi-tin) MAn) 



(23) 



with initial conditions Mq(0) — 1, M n> o(0) = 0. It is easy 
to integrate numerically this equation in order to have an 
idea of the solution. We find that the form 



Mi{n) = a 



(log(2n+l)) 1 



V2n + 1 

is solution, with = 3aj_i/(2(« — 1)), i.e. 



■> \ l 



(t-1)! 



(24) 



(25) 
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Fig. 4. Probability density function (PDF) of the occu- 
pation of each generation M after respectively (from left 
to right) 10 3 , 510 3 , 10 4 , 210 4 , 510 4 , 10 5 vortex splits. The 
PDF have been computed using the mean field equations 
(|23p . The shape appears as self-similar and is reminiscent 
of what is observed in Fig. [2J 
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Fig. 6. Probability density function (PDF) of the abso- 
lutes values of vortex strengths after respectively (from 
right to left) 10 3 , 5 10 3 , 10 4 , 2 10 4 , 5 10 4 , 10 5 vortex splits. 
The PDF have been computed using Ea. (|29l) . where we 
took into account only up to generation i = 40. The results 
appear to be coherent with what is displayed in Fig. [5] 
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Fig. 5. Maximum i* versus number of vortices N. i* is 
obtained from the distributions displayed in Fig. [4] A lin- 
ear scaling is observed as expected, with a measured slope 
of ~ 3.1, "close to the expected value 3/2 x In 10 « 3.45" . 

This can be checked easily by induction, we assume that 
Mj_i is of this form, and then we can solve Eq. (|2"3"]l for 
Mi in the continuous time limit in which it becomes 

dMj = M a _x _ Mj 
dn 2n + l 2n + 1 ' 

In this way we obtain that at a fixed level i, the occu- 
pancy Mi first increases with time, then decreases, with a 
maximum at n*(i) ~ exp(2i). At fixed time n on the other 
hand, Mi(n) has a maximum at i* ~ 31og(2n + l)/2. A 
numerical integration of the global populations given by 
Eq. (|2"3"]) is displayed in Fig. 0J One can notice similarities 
with the distributions of vorticity although the distribu- 
tions are more peaked. In order to test as well our anal- 
ysis, the location of the maximum of the distribution is 
displayed in Fig. [5] and a good agreement with the loga- 
rithmic law is found. 

Let us now compute the distribution of vorticity p{k, n) 
assuming we know the occupancy Mi(n). Hence let us con- 
sider a vortex living in the generation i. It has been the 
result of i splitting. Since the splitting rules (IT21 have no 



preferred order (as mentioned earlier they permute if we 
add 27r/3 to the random angle), we assume that the ob- 
tained vortex is always the third vortex hence its absolute 
vorticity k will end up being 

fc =4n( i - 2sin ^). ( 27 ) 

6 3 = 1 

and consequently its logarithm is 

i 

log =^log|l-2sin0j| -ilog3, (28) 
j=i 

with 9j being uniformly distributed random variables in 
[0 27r[. We can then gather the probability distribution of 
vortex strengths at generation i, which we note pi(k). And 
thus the vorticity distribution after N division writes 

p(k > n) = Y / Pi(k)p(M i> n), (29) 

i 

with p(Mi,n) = Mi/(2n + 1). 

In order to check these results, we compute p(k, n) 
using Eq. (f2T))) and compare the results to those displayed 
in Fig.f2J In fact, given the shape of the occupancy Mj(n), 
we can compute p(k, n) using only a "few" distributions 
Pi(k). For instance in Fig.[S]we computed the pi(k) taking 
into account in the tree the vortices only up to generation 
i — 40. We notice as well a very similar behaviour as the 
one displayed in Fig. [21 



5 Conclusion 

This paper is a first attempt at analysing the distribu- 
tion of vorticities originating from one point vortex using 
a dynamically compatible process, namely a reverse col- 
lapse route with the conservation of total vorticity. The 
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splitting process has been analysed in great details and 
shown to be dissipative. Afterwards a simple process con- 
sisting of randomly successive splittings is proposed and 
the resulting distributions have been analysed. Analyti- 
cal computation of the proposed process have been made, 
resulting for instance in the computation of the vortic- 
ity distribution after n consecutive splitting of vortices 
and show very good agreement with the numerical simu- 
lation of the process. This paper is a first step for further 
work. One could for instance modify the splitting rules 
in order to obtain a conservative process, but also could 
try to pick vortices and change how the splitting is done, 
with a non uniform probabilities, such as a Gibbsian one 
and analyse the resulting distributions. In other words, we 
could perform statistical physics of point vortices allow- 
ing a varying number of vortices according to the collapse 
rules, and see if this possibility changes the equilibrium 
features. Last but not least, it will be important to com- 
pare the obtained distribution with real data involving 
physical stochastic processes. In particular, it is tempting 
to compare the vortex distribution obtained with this pro- 
posed mechanism or its variants with experimental results 
obtained on two-dimensional physical flows. Remarkably, 
one has to notice that, starting with a positive vortex, 
the reverse collapse course induces naturally the creation 
of negative vortices, thus the engendered distributions will 
consist of both positive and negative vortices. Work is cur- 
rently under way to analyse these different possibilities. 
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